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Abstract: We discuss a new mechanism leading to a matrix product form for the stationary state 
of one-dimensional stochastic models. The corresponding algebra is quadratic and involves four 
different matrices. For the example of a coagulation-decoagulation model explicit four-dimensional 
representations are given and exact expressions for various physical quantities are recovered. We 
also find the general structure of n-point correlation functions at the phase transition. 
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Even for complicated one-dimensional many-particle models, the ground state can have a simple 
form. In spin problems it may be the tensor product of factors referring to single sites. While the 
correlations in this case are trivial, this is not so for a generalization where the product state is 
formed using matrices It has been found that the ground state of certain spin-one models 

and the stationary state for classical particles diffusing between two reservoirs have such a form. 
Also excited states have been described by the same Ansatz [^. So far, however, only diffusive 
systems could be treated successfully in this way. It is the aim of the present letter to show that 
the approach also works for more general situations. As an example, a particular reaction-diffusion 
model will be studied. 

We consider a stochastic two state model on a one-dimensional lattice with N sites. Its config- 
urations are defined by the occupation numbers ri , r2 , . . . , r^r each of which can take values and 
1 . We say the system has a matrix product ground state if its stationary probability distribution 
Po{ti,T2, . . . , ttv) can be written as 

N 

Po{ti,T2,...,tn) = ZN-UW\l[{TjD + il-rj)E)\V) (1) 

i=i 

where E and D are square matrices and {W\ and \V) are vectors acting in an auxiliary space. Z^v 
is a normalization constant defined as = {W\{D + E)^\V). The matrix product in Eq. (jl]) can 
be written formally as a tensor product so that the stationary state |Po) represented as a vector in 
configuration space is given by 



|Po) = Zr,'' {W\ (^g j \V) . 



The matrix product representation is a powerful tool since various physical quantities like the 
particle density 

(W\ C^^^DC^'^ W) 

faV= ' \w\c^\v) = " + (3) 

can be computed directly. Correlation functions are given by similar expressions in which C plays 
the role of a transfer matrix. 

The matrices used for the above Ansatz may by finite or infinite dimensional ^. We are 
going to study an example for the first case below. The fact that the probability distribution of 
some system is given by a product of finite dimensional matrices has far reaching consequences. 
Depending on the properties of the matrix C, correlation functions in such systems can have two 
forms which we want to discuss briefly at this point. Suppose first that the d-dimensional matrix C 
is diagonalizable and has eigenvalues Ai, with Ai < A2 < ... < \d-i < A^ . Then any n-point 
correlation function can be written as 



{Tj,Tj^...TjjN = J2 Civ({o-f }) exp { - ^ (^^.) ^ [ {jl - l)a'^ + (j2 - jl - 1)0-2 



+ (J3 - J2 - 1)4 + ... + (Jn - jn-1 " + (N - jnX+l ] } (4) 



where the first sum runs over all trf = 0, 1 under the restriction X]^=i < 1 . The quantities = 
{log(Arf/A^)}~^ are the correlation lengths. The CAr({(Ti}) are some coefficients which depend on the 
system size N and approach constant values for ^ 1 . All correlations depend exponentially on 
the distances involved |^] and the number of length scales equals the number of different eigenvalues 
of C minus one. 

The situation changes if C is not diagonalizable. In this case the matrix can be classified 
according to its Jordan normal form. As long as the Jordan block Jmax of the largest eigenvalue is 



1 



one-dimensional, the correlation functions again decay exponentially (the only difference to Eq. (^) 
is that algebraic prefactors to the exponentials may occur). On the other hand, if the dimension I 
of Jmax is larger than one, the correlations are dominated by algebraic terms with positive powers. 
One can show easily that the n-point correlation functions are given by 

{rnrj,-rjjN = E [^n{W^}) (ji - ir(i2 - Ji - l)'^^(i3 -J2 - 1)'^^ 

{(7i=0,...,/--l} 

■ (jn — jn-1 — 1)'^"(A^ — Jn)'^""'"^ ] + exponential terms . (5) 

Here the C7v({o"j}) are coefficients the large- asymptotics of which are generically proportional 
to A^^('^^) . The exponential terms are of type @ with algebraic prefactors and are generally 
negligible for large distances ji,(j2 — 3i), {jn — jn~i),{N — jn) ^ 1 • The correlations @ 
are completely different from those in a system with diagonalizable C. They involve the powers 
0, 1, 2, ^ — 1 of the distances. Only in special cases where all the C]sf{{ai}) vanish the correlations 
are of type (^). It is worth mentioning another special case: If / = 2 and the element Dd^d-i of the 
matrix D is zero, then all CN{{ai}) with more than one fij equal to 1 vanish, i.e. the correlation 
functions are linear in the positions ji,...,jn- If in addition -D^^d = 0, the algebraic part of 
the correlation function depends on the arguments (iV — jn) and only, i.e. {TjiTj2---'^jn 
c'(l — jn/N)+ exponential terms for N ^ 1 and (A^ — jn) ^ 1 where c' is some constant. Let us 
also note that the stationary correlations of any system with a ground state (|l]) containing finite 
dimensional matrices do not involve negative or non- integer powers of the distances. 

Up to now, matrix product ground states have been encountered in two situations. The first 
is found for models with the Hamiltonian H = in which the two-site interaction h itself 

already annihilates the ground state, i.e. hj j+i\0 >= 0. Here the algebra of the operators E and 
D is given by 



h 



0. (6) 



An example for this type of models is the spin-1 antiferrromagnet discussed in Ref. The second 
case is realized in models with open boundaries and particle input and output at the ends of the 
chain. These models are described by a time evolution operator H = hjj+i + h[^^ + h!'^^ 

where /i^^^ and h^^^ are 2x2 matrices for particle input and output. Here the basic mechanism of 
the matrix product ground state relies on the fact that application of /ij.j+i yields a divergence-like 
term on the right hand side 



h 



E\ ^ E 



(7) 



where e and d are numbers, normally e = —d = 1. Summing over the two-particle interactions, 
all these contributions cancel in the bulk of the chain. The remaining terms at the boundaries are 
canceled by a proper choice of the vectors (W\ and \V): 



< 



so that H\0) = 0. The most important two-state model of this type is the asymmetric exclusion 
process with external particle input and output |0]-[1,||6|. There are also three-state models to 
which the matrix product Ansatz has been applied |5[. But, as mentioned in the beginning, all 
known examples are diffusive systems. 

The generalization which we are going to use, consists in replacing the numbers e and d by 
matrices E and D. The idea goes back to Ref. where the special case E + D = has been 
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introduced in order to solve the time-evolution of the asymmetric diffusion model in one dimension. 
The generalized algebra 



< 



h 



W\ h^^^ 







< W\ 







W >-- 



\V > 



(9) 



(10) 



is quadratic on both the left and the right hand side. In contrast to the usual matrix Ansatz (|^) 
the generalized Ansatz (|9|) can be applied to systems which include particle reactions. 

As an example we consider the asymmetric coagulation-decoagulation model. In this model 
particles diffuse on a linear chain. When two of them meet, they can merge (coagulate) to a single 
one. In the same way a single particle can split up (decoagulate) into two particles. Assuming no 
particle input and output, we therefore have six different processes: 



diffusion to the left 
diffusion to the right 
coagulation at the left 
coagulation at the right 
decoagulation to the left 
decoagulation to the right 



+ 

^ + 0- 
A + A 
A + A 
iD + A- 
A + 9- 



^^ + 
^0 + ^ 
^^ + 
^ + A 
■A + A 
■A + A 



with rate ai 

with rate aji 

with rate cl 

with rate cr 

with rate cIl 
with rate 



In what follows we consider the special choice ai = cl = q, an = cr = q~^, = Ag and 
dR = Ag^^ where the diffusion and coagulation rates coincide and all reactions have the same bias 
in one spatial direction. Since in this case the model can be mapped on a free fermion model, it 
is integrable and various exact results have been obtained [0, ^- The model is controlled by 
two parameters, namely the asymmetry parameter q and the effective decoagulation rate A. Its 
phase diagram shows two phases, a low-density phase for A < — 1 and a high-density phase for 
A > — 1. At the phase transition point A = — 1, the gap in the relaxational spectrum vanishes 
and algebraic long-range correlations can be observed In a basis (00, 0^, ^0, AA) the two-site 
term in the time evolution operator H 



Z^7 = l 










(A + l)g 

-q 

-Aq 



hjj+i reads 




(A + l)g-i 
-Ag-i 





-q' 



q + q'^ J 

EE - 
ED 
DE 
DD 



Therefore the bulk algebra (P) is given by 



{A+l)qED - q-^DE - 
-qED + {A + l)q'^DE 



AqED - Aq'^DE + (q + q 





qDD 
)DD 



EE 
ED 
DE 
DD 



and the boundary conditions (|lO|) read 

{W\ E = {W\ D = E\V) = D\V) = 0. 
Writing C = E + D, (7 = ^ + L) and 7^ = A -M, the algebra (|l2|)-(|l5|) simplifies to 

[C,C] = [E,E] = 



EC 
CE 



EC 
CE 



{^\ + q~^)EC 
{7\-'+q)CE 



YqEE - q-^CC 
-f^q-^EE - qCC. 



(11) 



(12) 
(13) 
(14) 
(15) 

(16) 

(17) 
(18) 
(19) 
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In contrast to algebras for diffusive systems (0)-(|8|), tfie above commutation relations do not allow 
to reduce the number of factors in a given product of matrices. Therefore products of different 
lengths are independent. Products of the same length, which correspond to a given system size, 
obey linear relations as follows. For a given product, e.g. CECEE, we compute the difference 
CECEE - CECEE by using the commutation relations (0)-(|ll). Writing {W\... \V) = (. . .) 
and using {CECEE) = {CECEE) = one obtains: 

-1 



{CECEE) 



In general, if {Pi^\ ^ € 1 . . . N]^} is the subset of products with factors containing k matrices C, 
linear relations of this type have the form 

Nk+i 



+ 



■q-^ {EEC EE) 



(20) 



+ 7^iq + q~^){CEEEE) + {q + q-^) {CCCEE) + q{CECCE) 



+ E 



k,k-l I Tj(k-l) 



l...iV- 1) 



(21) 



As can be seen from the commutation relations, the coefficients c^'j^^ obey 

Nk±i 



< f < 1 



n ^ k,k±l ^ 



k,k±m 



(22) 



E E<r" = i 

j=l m=±l j=l 

Therefore by iterating Eq. (21) one gets more and more complicated linear expressions with positive 
coefficients which involve all subsets k = . . . N. Since there are no such relations for A; = and 
k = N, one finally ends up with only two contributions: 



{E^) + {l-af^){C^). 



■<N\ 



(0 < af ) < 1) 



(23) 



The expectation values {E^) and {C^) are independent. Therefore the vector space of words 
of a given length decomposes into two subspaces in which the expectation values are proportional 
to {E'^) or (C^) , respectively. Consequently physical observables are parameterized by the ratio 
A := {E^) / {C^) . This is related to the fact that the model has two independent ground states, 
a trivial one which is the empty lattice (A = 1) and a nontrivial one where particles are present 
(A = 0). 

A trivial representation of the above algebra (p!6[)-(p!9|) isi? = C= l, £' = (7 = which describes 
a system without particles. In the symmetric case q = 1 there also is a second one-dimensional 
representation E = 1, C = 7^, E = C = corresponding to a factorized ground state with finite 
particle density A/(l + A). In the general case q ^ 1 the model is known to involve three different 
length scales, and therefore any nontrivial representation of the algebra has a dimension d > 4. 
Furthermore representations of the algebra may be different in each sector so that they may depend 
explicitly on and A. We found a four-dimensional representation which is given by 



El 
















1 





7" 






^2 



El 






Vo 



Aiq-q 




1 

-1^ 



-Q 
-Aq 

/ 



Ci 







( q-' 


q~^ 






Ci 







1 1 












7^ 















q^ ) 








Aq-^ 


g-1 




1 _ 





A{q 


-'-1) 


q-q ^ 






















vo 
















(24) 
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{Wi\ = (l - g2, 1, 0, a) , \Vi) = {b, 0, - 1 



where 



-2Ni 



{q' - l') + (7)-'^(7' - + 1) - q-''"" {q'l' - 1) 



-2N (^2^,2 



(^2 _ q2)(-^2 



-2) (^2 



(25) 



+ 



A- 1 



(7^ -9^) (7' 



2Uo/2 _ ^-2W^2 



)(g2-^-2) 



The case IS. ^ q^ — 1 : For practical purposes it is desirable to have a representation in which the 
matrix C is diagonal. For IS. ^ q'^ — 1 an appropriate similarity transformation yields 



/ q-^ q^-T^ q^-i '7'(l-7') \ 



E2 

{W2\ 
\V2 



I 



,-2 






V 



(I 



1 



gi2-j,2 



7 




, 0. 





1 





T -q 

l\q^ - 1) 



a{q 



Co 



J 



\ 



\ 

-'){l'-q'h' -q'l'' 



) - ( 



g2^2 _ I ' 

6((?4-l)(g2^2_i) + ^4 



(72-l)(g2 + l) 

?2(^2 _ 1) (^2 _ 1)^2 



(26) 



— 'j^q'^ 



g2 + 1 'j^ 

Using this representation, it is easy to derive the particle density (^) in the sector A = 

72^((7^ - 1) + iq' - ihHq^)-'') - ^'^((7^ - i)q'-'' + iq' - i)iqhr'') 



N 



^2 (^2Af _|_ ^- 



■2N _ g2N _ q-2Af j 



(27) 



which coincides with the result obtained in Ref. [Q]. We also checked that the two-point correlation 
function {TiTj)^ is obtained correctly. 

The case A = q'^ — 1 (q > 1) : Here the two largest eigenvalues of the matrix C, namely 1 + A 
and q'^, coincide and C is not diagonalizable. We therefore choose a representation where C has 
Jordan normal form 



f q ^ q^ — q 
g-2 

V 



2 _i g^+g-2-2g2 \ 




1 




q-' 

l+2q'^ 



( 



\ 



1 
g2 

\ 





1 

^2 



(28) 



{W^\ = [q- 
where 



1, 0, i + 



\V^) = {l + h{q^ -q-^), 0, -1, q^-q- 



2N -2Ni 
q a — q 



LA(g2L_^-2L) 



Lq 



-2L 



-,2L 



-2L 



1-A - r - 

Using this representation, the density at a site j in the sector A = is easily obtained as 



-4j 



(29) 



(30) 



which agrees with the result from Ref. There is no term proportional to (j — 1) or (j — 1){N — j) 
because of ,04,3 = -D4,4 = (see Eq. ([2^)). It turns out that any n-point correlation function 
{'^ji'^j2---'^jn)N depends only on the two positions ji and jn ■ According to our discussion at the 
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beginning of this letter, its algebraic part is a linear function in j„ only. In fact it is given by 
^OTZT (1 - 9"^)"""^ {(1 - (l - + ^} • The exponential part of {Tj^Tj^...TjjN decays 
with ji on length scales {2logq)~^ , (41ogg)~^ and with j„ on the length scale (2 log (7)^^ . 

While an Ansatz of type (|^) with an algebra (^), ([l^) can be made for any one-dimensional 
reaction-diffusion model, it is not clear under which conditions a matrix representation really exists. 
In particular, we do not know if the existence of representations is related to the integrability of 
the system. One should therefore investigate non-integrable examples. Also the extension to 
systems with open boundaries would be of interest. However, since some open systems are known 
to have correlations decaying with negative powers of the positions, the corresponding matrix 
representations are expected to be infinite-dimensional. 
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